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Abstract 

We present new analytic solutions to the relativistic Boltzmann equation within the relaxation 
time approximation. We first obtain spherically expanding solutions which are the kinetic coun¬ 
terparts of the exact solutions of the Israel-Stewart equation in the literature. This allows us to 
compare the solutions of the kinetic and hydrodynamic equations at an analytical level. We then 
derive a novel boost-invariant solution of the Boltzmann equation which has an unconventional 
dependence on the proper time. The existence of such a solution is also suggested in second order 
hydrodynamics and fluid-gravity correspondence. 

TAGS numbers: 25.75.-q, 12.38.Mh, 24.10.Nz, 51.10.-ty 


1 


I. INTRODUCTION 


The Boltzmann equation is one of the most important equations in contemporary physics 
that allows us to characterize the transport prOTerties of a dilute gas based on the micro¬ 
scopic dynamics of the constituent particles [l, 2|. It is a nonlinear partial differential 
equation for the one-particle distribution function f{x,p) which is very difficult to solve ex¬ 
actly by analytical means. Due to this limitation, it has been more convenient to study the 
equation by solving it either numerically or by developing approximate solutions based on 
different expansion schemes. Both approaches have been successful when comparing their 
predictions with experimental results, but at the same time, they have their own limitations 
especially in a relativistic setting. Thus, exact solutions to the Boltzmann equation would 
certainly be useful to constrain the validity of different perturbative and numerical methods. 

Symmetries provide powerful methods to solve and simplify complex problems in physics. 
This has been particularly useful in relativistic kinetic theory. For instance, for the rela¬ 
tivistic Boltzmann equation within the relaxation time approximation (RTA) j^. Refs, 
considered a solution having the same symmetry as the Gubser flow 6| which is a boost- 
invariant solution of the conformal hydrodynamic equations relevant to heavy-ion collisions. 
The symmetry group of the Gubser flow restricts the number of independent variables as well 
as their particular combinations on which the distribution function depends. As a result, the 
Boltzmann equation can be effectively reduced to a one-dimensional problem and an exact 
solution has been constructed j^, By using similar symmetry arguments, other solutions 
of the RTA Boltzmann equation have been found in the literature for near equilibrium BJii 


and highly anisotropic systems 


14 


17|. In these results, the solutions are written formally 


in terms of the effective temperature of the system which has to be determined numerically. 
While this can be done straightforwardly in practice, it is always welcome to have fully an¬ 
alytical solutions where one can understand various aspects of nonequilibrium dynamics in 
a completely controllable manner. In this paper, we make progress in this direction by pre¬ 
senting new analytical solutions to the RTA Boltzmann equation for conformally invariant 
systems. Each of the solutions is shown to have an explicit counterpart in relativistic viscous 


hydrodynamics discussed in 


18 


19l |. This allows us to not only compare the solutions of 


kinetic and hydrodynamic equations at an analytical level, but also shed light on how the 
hydrodynamic solutions are obtained as a coarse-grained version of the kinetic solutions. 


2 























This paper is organized as follows: In Sect. |TT]we describe the kinetic theory approach to 
the Hubble flow solution Q3- In Sect, [ml we derive a new boost-invariant solution to the 


RTA Boltzmann equation which features an unusual dependence on the proper 


ime. We 


provide evidence that such a solution can exist in second order hydrodynamics and also 
in fluid-gravity correspondence 2l|]. The conclusions of this work are presented in Sect. IIVI 


II. KINETIC THEORY DESCRIPTION OF THE HUBBLE FLOW 


Our starting point is a spherically expanding solution of the relativistic ideal hydrody¬ 
namic equations which is characterized by the following flow velocity and the energy 


density e 


20 | 


= 


t r 

Tr ' Tr 


e oc 


T 


( 1 ) 


where Tr = — 'P (r = {xi, X 2 , x^)) is the proper time. The solution ([T]) is valid for 

the conformal equation of state e = 3p {p is the pressure) which we assume throughout this 
paper. It is convenient to switch from the Minkowski coordinates to the following coordinate 
system via a Weyl rescaling 


ds'^ = —dt^ + dr"^ -I- r‘^{d6‘^ + sin^ Odcj)'^) 

ds^ = —5- = —dx^ + dv^ + sinh^ v{d6‘^ + sin^ 9d(j )'^), 


( 2 ) 


where y = Inr^ and we introduced the ‘rapidity’ variable v = tanh“^ j. In this coordinate 
system, the flow velocity is simply = 5^^ and the energy density i = r^s is a constant. 
We shall refer to this solution as the three-dimensional (3D) ‘Hubble’ flow in analogy to the 
well-known flow solution in cosmology. 

Our goal in this section is to describe the 3D Hubble flow and its nonequilibrium gen¬ 
eralizations within relativistic kinetic theory. The relativistic Boltzmann equation for the 
distribution function of massless particles f{x,p) in any curved spacetime reads as QQ 


+ r^iPAP" (/ - fe,) , 

^ OPi Trr 


( 3 ) 


where we employed the so-called relaxation time approximation (RTA) [3[ in which the 
collision term is linearized around the equilibrium distribution feq to be specified shortly, r.^- 
is a characteristic time of the order of the time between succesive collisions and in general, 
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it can depend on the space-time and momentum-space coordinates. In the above equation, 
/ = f{x^,Pi) is considered to be a function of the space-time coordinates and the three 
dimensional spatial momentum components Pi {i = 1,2,3) with lower (covariant) indices^. 
The energy of the particle is determined from the on-shell condition g^^p^Pu = 0. 

We analyze the equation ([2D in the coordinate system x^ = (y,r;,0,0). In this case. 
Pi = {pv,Pe,P 4 >), and the on-shell condition becomes 


2 2 I 

Py=Pv + 


Pn 


' • 1,2 ’ 
smh V 

^2 — ^2 , ^2 / „;„2 


(4) 


where we abbreviated Pq=pI+ sin^ 9. Computing the Christoffel symbols and using 
the flow velocity it^ = 6^, we hnd 

(p^d^ + p^d^ + p^dg + p'^d^ + f = -^P^if - feq) ■ (5) 

V smh'^ V opv smh^ v sm'^ 9 opg J 

where we dehned (5 = TrjTT,? In the following, we assume that / is independent of 0. 

First let us specify the equilibrium distribution. Knowing that the flow is static in this 
coordinate system, we immediately hnd that the Boltzmann distribution 

/ = /„(?>=) = , (6) 

^ Since / is 


where the temperature T is a constant and pP^ is as in (|4D, exactly satishes 
a scalar invariant, in Minkowski space the equilibrium distribution is 


feq 6 


(7) 


where = p^/xr and T{Tr) = T/xr is the temperature in Minkowski space. Thus, the 
distribution function (ED is the kinetic counterpart of the ideal Hubble how. 

We now add perturbations on top of the ideal solution. In a conformal theory, oc 1/T 
by dimensional analysis. This means that oc r^, and therefore /? in Eq. ([5D is a constant. 
Writing f = feq + Sf = feq{^ + *h), we hnd the following equation for <F 


p^dy^ p^'d^ + p^dg + 


Pq cosh V d pjj, cos 9 d 


+ 


sinh^ V dpv sinh^ v sin^ 9 dpg 


$ = -/3pX$. 


( 8 ) 


^ The Boltzmann equation ([Sj is covariant in coordinate space but not manifestly covariant in momentum- 
space [^ . One can sort out this problem by considering a manifestly covariant Boltzmann equation for 
off-shell distribution functions /(a;^,p^) [^. We will not consider this approach in our work. 

^ Under the Weyl transformation ([2]), Ttt is rescaled by a factor of [3. 

^ We may also take the Bose-Einstein or the Fermi-Dirac distribution as the equilibrium distribution. 
Actually, in the RTA any function feq{p^) satisfies ([5|). We need the full Boltzmann equation to uniquely 
determine the equilibrium distribution. 
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We shall solve (jH]) in two interesting cases corresponding to the scalar and tensor perturba¬ 
tions around the ideal solution. 


A. Scalar perturbation 


Motivated by a recent work 


24| . let us hrst discuss the solution which is independent of 


‘time’ x- Assuming 0(3)-symmetry, we hnd the following exact solution of 

4. 


2/5 — p„ tanh v 

= /C(p^,pQ)exp ( -itanh"^ 


( 9 ) 


The only constraint for the function K, is that it must satisfy the Landau matching condition 
for the energy density e = s/t^. 


£ = 


[u-p] /, 


( 10 ) 


(27r)3 y (27r)3 y 

where d?p = dp^dpedp^. If we further assume that /C does not depend explicitly on then 
the condition reduces to 


dp^p^e = 0 . 


( 11 ) 


The above solution (I9l) has been derived in an analogous way to 2^ based on a different ideal 


hydrodynamic solution. As discussed in that reference, this type of solutions is characterized 
by certain scalar moments of / and the associated entropy production despite the vanishing 
shear-stress tensor 

1 f d^p 


= 


= 0 . 


( 12 ) 


(2!r)3J 

where + AJA“) - jA'-'-A^^ with A'"' = g)^ + ti"!!". (E) means that the 

solution does not allow for a hydrodynamic description. Rather, it represents the relaxation 
of ‘fast’ degrees of freedom usually not taken into account in hydrodynamics. 


B. Tensor perturbation 


We now return to (0) and consider the tensor perturbations. We parameterize the 
nonequilibrium part in such a way that the connection to hydrodynamics is transparent 


2 2 
$ = V, 9) = ^ 


8 T 6 


8 T 6 




(13) 
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The first equality is the standard parameterization in the moment method where is 
identified with the shear stress tensor in viscous hydrodynamics flT^ . In the second equality, 
we extracted the exponential relaxation factor which turns into a power-law behavior 
in the original Minkowski space. 

Let us hrst look for 0(3)-symmetric solutions where $ is independent of 6 and = 
Substituting Eq. (IT^ in Eq. (|8]), we hnd the following equation for 


P^'d^+Pn- 


cosh V d 


sinh'^ V 


pI 


pI 


2 sinh^ V 


= 0 - 


(14) 


However, we immediately encounter a difficulty. It is easy to see that there can be no 
solution to flTT|) . Indeed, the general solution of the differential equation 


P'^dv+Pn- 


cosh V d 


$ = 0 , 


(15) 


' sinh^ V dp^ 

is $ = ^{p^,Pq), and this is clearly incompatible with the structure of flTT]) . 

The trouble is that this negative result is in apparent contradiction to the finding in 


Ref. [19[. There, the authors derived exact solutions to the Israel-Stewart equation in 
hydrodynamics which relax to the Hubble flow at large times. In the present notation, 
the 0(3)-symmetric solution is (see (74) of [l9j) 

1 


ttL, oc 


(16) 


sinh^ V ’ 

which, however, is not a solution to (fTTj) as we have just seen in the above derivation. On 
general grounds, it is expected that every solution of the hydrodynamic equations has a 
microscopic counterpart in kinetic theory, and actually our motivation here is to rederive 
(1T6|) as the solution of the Boltzmann equation. It is tempting to think that this is a problem 
of the RTA which oversimplihes the collision term of the Boltzmann equation and therefore 
restricts the solution space of kinetic theory. Yet, one can derive the Israel-Stewart equation 
starting from the Boltzmann equation in the RTA (see for example 25(]), and this implies 
that the above conflict must somehow be reconciled within the RTA. 

In order to understand in what sense Eq. (1T6|) is a solution, let us substitute it into 
Eq. dUl) 


Pvdv + Pq 


cosh V d 
sinh^ V dpi 


pI 


1 —pv cosh V 

2 sinh^ V J sinh"^ v 2 sinh"^"*"^ v V 


Pn 


(2Ap: 




sinh V 
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where A = ?>. As expected, the right-hand-side does not vanish for any value of A. However, 
the question is whether these unwanted terms affect the hydrodynamic equations. The 
energy-momentum conservation equation is 

1 


0 = = 


(2i)s 




(18) 


Taking the component = v, we see that all that is needed for the hydrodynamic equation 
to hold is that the following integral vanishes 


d^p 


:Pv (Eq.([IZ])) oc 


d^p 


pI ^ApI 


{A + 6)pI 


cx A — 3 


(19) 


J V sinh^y 

This is indeed the case when A = 3. It is easy to see that the components of the equation 
flTSll other than p = v are trivially satisfied even when A 7 ^ 3. 

The above analysis teaches an important lesson about comparing solutions of kinetic 
and hydrodynamic equations. Since the hydrodynamic equation is a course-grained version 
of the kinetic equation, it admits a class of solutions which are not sensitive to the exact 
details of kinetic theory. To accommodate this, we have enlarged the solution space of the 
Boltzmann equation in the RTA to allow the equation to be satisfied up to terms that do 
not affect the macroscopic (hydrodynamic) equations. With this qualihcation, remarkably 
the solutions of the Boltzmann and hydrodynamic equations are exactly the same in the 
sense that they are characterized by the same macroscopic variables e, tt^^, etc. This is in 
contrast to the common perception that solutions of the hydrodynamic equation are only 


an approximate version o 


in the literature 


4 i5i, m 


the so’ 


15 


25 


utions of the Boltzmann equation, as repeatedly observed 


26|. 


In Ref. [19|], along with the 0(3)-symmetric solution (ITbl) . non-0(3)-symmetric solutions 
to the Israel-Stewart equation were also obtained. It is straightforward to generalize the 
present analysis to this case. Here we consider only one of the non-0(3)-invariant solutions 
found in which reads'^ 


^ -^ 

V ~ 9 ~ 4 ) ~ sT ^ 173 ’ 

2 ^ sm U smh v 


( 20 ) 


We have checked that the other non-0(3)-symmetric solution in 19|, which is proportional to ^^^ 3/2 ^ ; 

also satisfies the RTA Boltzmann equation following exactly the same pattern. 
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so that 


(v 6) oc. { -\ _ — _^ ^ _ 

^ sinh^ V sinh^ v sin^ 9) sinh^ v sin^ 9 


( 21 ) 


Substituting into ([8]), we find 


Pocoshn d p‘icos9 d . „ 

P d^+ p^de + ^ , 3 —tt— + —ri— . 3 ^ ) P P 


oc 


sinh^ V dpy sinh^ v sin^ 9 dpe 


pecot9\ I 2 , 
Py cosh V H-—- Pu + 


pI 


^pI 


( 22 ) 


sin^ 9 sinh"^ v v J sinh^ v sinh^ v sin^ 9 , 

As it happened in the 0(3)-symmetric case, the right hand side of (|22|) does not vanish. 
However, now we know how to sort out this. The nonvanishing terms in Eq. fl22|) do not 
affect the hydrodynamic equation (flSD because 


cPp 


:Pv (Eq.([22])) = 


d^p 


-pe (Eq.([22])) = 0, 


(23) 


J V^p^ ' J 

for p = v,9. (The other components are trivial.) In this sense, Eq. fl?I|) is the solution to 
the Boltzmann equation which corresponds to the hydrodynamic solution (l20ll . 

Finally, we note that in both the 0(3)-invariant and non-invariant cases one can obtain 
the free streaming solutions by taking the limit /3 —?■ 0 

,2 


TT 


/ = /.q 1 + 


(24) 


In the absence of the exponential damping, the system never reaches thermal equilibrium 


(see, also 


. t 


23) 


III. BJORKEN FLOW REVISITED 


In the previous section, we presented explicit analytic solutions of the Boltzmann equa¬ 
tions in the RTA. This has been possible largely due to the fact that the ideal hydrody¬ 
namic solution is static in a cleverly chosen coordinate system. For essentially nonstatic 
flows, analytic solutions of the Boltzmann equation are difficult to obtain even in the RTA. 


Nevertheless, in the boost-invariant case relevant to heavy-ion collisions 


, one can gam 


analytic insights into the behavior of the flow at late times. In this section we revisit this 
problem and elucidate a new solution of the Boltzmann equation which exists in the presence 
of conformal symmetry. 

















We work in the coordinate system 


ds^ = —dr"^ + + d^T + x^dcf)"^ , 


(26) 


where xt = ^Jx\ + X 2 is the transverse coordinate, r = and C, = tanh“^ ^ are 

the one-dimensional analogs of the three-dimensional proper time and the rapidity v 
introduced in the previous section. The RTA Boltzmann equation with the comoving flow 
velocity = 5^ 281 was first studied in j?] and developed more recently in js-lS]. Assuming 
that / depends only on r, we need to solve 


drf = -(/-/e,) 

'^TT 


(26) 


The equilibrium distribution is taken to be the Boltzmann distribution as before feq = 
^u-pjT _ g-p^/T(T) pT — ^p 2 unlike the Hubble flow case, T(r) is an 

unknown function which is dynamically determined from the first moment of the Boltzmann 


equation fl2^ which gives us the dynamical Landau matching condition 


Q 


e = 


( ^2 f 1 

[u ■ p) feq = 


dfp 


{u-pf f. 


(27) 


(27r)3 y (2vr)3 y yf^p'^ 

Consequently, / = feq does not solve fl26l) exactly, which is a manifestation of the nonstatic 
nature of the geometry. 

In terms of the energy density fl271) . the solution of fl26l) is formally given by QQi 


fjq-i 1 / q -'2 arctan 

e(r) = D(r, ro)£o(r) + j 7^ + “ 


- 1 


(28) 


- 1 


where 


D{t,to) = exp - 


'TO 


dr' 

T,(t') 


(29) 


and 


^o(r) = 


d?P, r^2 


( 2 vr) 


rp'’ 


{p") f{ro,Pc,PT)- 


(30) 


In Eq. fl30|) . /(tq) is the initial distribution at Tq , whereas p'^ = + Pt (on-shell 

condition) is defined at r. After integration by parts, Eq. (j28]) can be written as 




dTD{T,T')-^ (^e{T')R = D{t,To) (^eoir) - e(To)R , (31) 
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where we abbreviated 


R{x) 




(32) 


Following 
(1^ becomes 


0 . 


let us first assume that is a constant. Then = e 




, and 


^ (^e{T')R (£o(r) - e{To)R . (33) 

Since the right-hand side decreases with time as (9(l/r), one should not allow the left-hand 
side to grow exponentially in r. This leads to the condition 


A 

dr' 



= 0 . 


r'=T 


(34) 


With R{1) = 1 and R'{1) = 4/3, fl3T|) gives e(r) ~ 1/r^/^. Thus the Bjorken solution is 
recovered 28(]. The same conclusion is reached if Ttt depends on time as T 7 r(r) ~ with 
0 < p < 1. However, when p = 1, this argument breaks down.® In fact, there exists a novel 
asymptotic solution of the form 


e(r) 



(,r > To) 


(35) 


whose normalization constant C is not arbitrary but is an intrinsic parameter of a given 
theory. 

In order to verify this statement, we first note that the two conditions p = 1 and e 1/r^ 
are naturally related in the presence of conformal symmetry. Indeed, in a conformal theory, 
Ttt oc by dimensional analysis. We can then write in the asymptotic regime the following 


so that p = 1 together with 


Hr 


fn Tr_ T 

jl/4 “ - p ’ 


D(t, r') 



(36) 


(37) 


® The case p = 1 was previously studied in [sj without assuming conformal symmetry. It was found that 
the Bjorken solution e 1/r^/^ is recovered only in the limit ^ 00 (see Eq. (IMF , whereas for other 
values of /3 the asymptotic behavior is e 1/r'’' with 1 < 7 < |. The special solution (l35l) which we are 
going to derive was not noticed in [^. 
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lead us to write Eq. (^T\i as 


'TO 


( 7 )) “ “ e(r„)R 0 


(38) 


We see that, instead of an exponentially growing factor as in flH^ . the integrand contains 
only powers of r'. At large times r tq, the right hand side is of order C>(l/r) + C>(l/r^) 
with To-dependent coefficients. This should match the contribution from the lower bound 
t' = tq of the r'-integration on the left hand side. Then the contribution from the upper 
bound r' = r of order which is independent of tq , must vanish. 

For generic values of P, the r'-integral in Eq. fl55]) cannot be done analytically. In order 
to study the behavior near the upper limit, we expand R{t'/t) in powers of r — r' 

2 

+ ■ ■ ■ , (39) 


and integrate over r' term by term using e(r') ~ In fact, we need to expand to all 

orders in r — t'. In practice, we used Mathematica and expanded R{t'/ r) to 0{{t — 

We then require that the contribution from the upper bound r' = t vanishes. This yields 
the value® 


(3 ^ 1.27672. 


(40) 


In Appendix we present another derivation of this constant. We thus conclude that, 
in the presence of conformal symmetry, the Boltzmann equation in the RTA admits the 
following asymptotic solution 


C 


e = 


with ^ 1.28f^ 




(41) 


As announced, the normalization is completely determined by which is an 

intrinsic parameter of a given theory. Note that in QCD at high temperature which is 
nearly conformal, we have Tjt r]le (TA^lni) ^ (A = is the’t Hooft coupling) [29| 
so that 


s ~ 




(42) 


® In addition to this, Mathematica finds other positive roots such as /3 = 2 and (3 = A. We discard them as 
artifacts. For these values of /3i one can evaluate the r'-integral explicitly, and find inconsistencies with 

dMH- 
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We now discuss the connection to hydrodynamics. In fact, a solution very similar to 
Eq. fHTj) has been previously found in Ref. as an exact solution of the general second- 
order hydrodynamic equations and dubbed ‘unorthodox Bjorken flow’. This is given by (see 
Eq. (112) of Ref. (idI) 


E = ^ with = 1th - S’) - 2^. 




6 


(43) 


where rj = is the shear viscosity and = e is one of the second-order transport 

coefficients in the constitutive equation for 


TT' 


/if/ _ 




_Ol\Q 


(44) 


1/3 V r ” ' 7r7r^Q!/3'' 

These transport coefficients can be evaluated from the Boltzmann equation in the RTA.^ 
Using Ttt = 57]/(Ts) = ISr//(4^) from the Chapman-Enskog theory (s is the entropy density) 


25 


and Tt^tt = 10r7r/7 for the massless Boltzmann gas 

= f^TT, 2.06 f, 


30|, we get 


(45) 


which is in the same ballpark as Eq. fl4T|) . The discrepancy may be alleviated in a more 
precise evaluation of these coefficients from the Boltzmann equation. Conversely, if we 
assume that the Chapman-Enskog value for rj is precise, our result may be used to estimate 
the value of in the RTA. Equating (l4Tll and (03]), we obtain r,r 7 r ~ 3.77 Tt, , which 


out to be a few times larger than the results previously derived by different authors 


25 


urns 


30|. 


A. Fluid-Gravity duality 


In this subsection, we point out that the existence of the solution of the type (BBi 
is also suggested in the framework of fluid-gravity correspondence.® In strongly coupled 
A7 = 4 supersymmetric Yang-Mills theory, hydrodynamic flows are dual to solutions of the 
five-dimensional (5D) Einstein equation in asymptotically Anti-de Sitter (AdS) spaces (see 


Ref. 


21 | for a pedagogical review) 

Rab 


-GabR 


QGab — 0 , 


(46) 


^ The other second-order transport coefficients which are not derivable from the Boltzmann equation in the 
RTA have been ignored in (|45)) . 

® This subsection is largely motivated by interesting discussions with Jorge Noronha to whom we are 
grateful. 
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where Rab and R are the Ricci tensor and the scalar curvatnre, respectively. In the 
Fefferman-Graham coordinates, the 5D metric Gab can be written as 

dsl = GABdx^dx^ = ^ \g^u{x^ z)dx^dx'' + dz"^] , (47) 

where 2 ; is the ‘fifth’ dimension. Near the Minkowski boundary z = 0, the solution of fl46|) 
takes the form 

= gf} + z^g]$ + z^gf} + • • •, (48) 

where is the four dimensional fiat Minkowski metric and is related to the energy 
momentum tensor of in the gauge theory as 


T = 

27r2 


/V2 

(4) 
y flu 


(49) 


For boost-invariant flows, the dual geometry at asymptotically large r turns out to be 
the scaling solution [31 1 

1 


ds^ = — + dz^) 


(50) 


The scaling variable is defined by u = zjRl^ where s is a parameter required to satisfy 
0 < s < 4 from the positivity constraint of the energy density. The Bjorken solution 
e ~ corresponds to s = 4/3 and 


aiu) 


(1 - fuy 


b{u) 'zz c{u) 


1 I 4 

l + -« 


(51) 


(l + fn^) ’ 

where Cq is a dimensionful normalization constant. flHT]) satisfies the Einstein equation up to 
terms subleading (at fixed u) in inverse powers of r which are related to the viscous effects 


32 


. It was shown in 


3l| that within the range 0 < s < 4, the value s = 4/3 is the only 


acceptable choice after requiring that there are no singularities in the bulk. The boundary 
value s = 4 was not studied in |m| because, as the authors pointed out, their method to 


solve the Einstein equation (expansion in negative powers of r at fixed u) fails when s = 4. 
What happens at s = 4 is that u = z/t is dimensionless, and therefore at fixed u negative 
powers of r cannot appear. In other words, one has to solve the Einstein equation exactly. 
Furthermore, it is interesting to note that not only eo becomes dimensionless when s = 4, 
but also it becomes a fixed number due to nonlinearity of the Einstein equation as shown 
below. 
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The reason why we are interested in the special value s = 4 is that, if a solution with 
s = 4 indeed exists, then it is dual to a hydrodynamic flow e ^ ^ which has the same 
r-dependence as gHi. Unfortunately, we have not been able to solve the Einstein equation 
for s = 4 exactly in a closed form. Instead, we performed a perturbative expansion of a{u), 
etc. in Eq. fl50ll to high orders. The result is 


-u 


9 4 

f _ ^8 _ 4^10 

3 3 


_ 


-u 


^ ^ ” 33 2136 21 630 

^ o 4 10 fj ^ c 22 12 95 14 287 m 880 ^ 15761 

b{u) = 3u^ + —u^ + 3 m « + 4 ^ 1 ° + 

/ x 4 4 11 o 10 in 17 12 208 14 635 m 223 ^ 36661 sn 

c(u) = -u^ - -u^ - 2 m®----- 

^ ^ 3 3 3 21 36 7 630 


+ •••, (52) 


m 2 ° + 


H-. 


In fact, the expansion can be carried out up to all orders SJ]- Remarkably, to lowest order 
the coefficient of in a{u) is arbitrary, but once we go to higher orders the nonlinearity 
of the Einstein equation selects the coefficient to be —1. Assuming that fl52l) represents the 


near-boundary behavior of a well-defined solution, we read off the energy density in field 
theory 


e = 


2'k'^ 


(53) 


Once again, we hnd that the coefficient of l/r"^ is uniquely determined, this time due to the 
nonlinearity of the Einstein equation. Moreover, the energy momentum tensor has the form 

(64) 


c = = _rpxx _ _j^yy _ yCC 

3 


which is exactly the same as that for the unorthodox Bjorken flow 
comparison clearer, we note that by using Eq. 


19]. To make the 




We then use the known results at strong coupling [35|, [36 1 


£ = 




Tn = 


2 - ln2 
27iT 


to rewrite Eq. 


as 


e = 


{f3f^ 


with = 


( 55 ) 


(56) 


(57) 


T' 2 - In 2 

which is rather close to Eq. gB, though of course the agreement is not expected since one 
is comparing strongly coupled and weakly coupled theories. 
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Unfortunately, we have not been able to resum the series (l52|) in a closed form and this 
prohibits us from studying the possible singularities of the geometry in the bulk. Moreover, 
the radius of convergence for the series 0521) seems to be decreasing towards zero. Therefore, 
at the moment we cannot draw a conclusion about the existence of solutions of the form 
although it would be very interesting to have one in view of the results in kinetic theory 
(ITT]) and hydrodynamics (H3|) . We leave this problem to future work. 


IV. CONCLUSIONS 

In this paper, we have derived analytic solutions of the Boltzmann equation in the RTA for 
conformally invariant systems. In the ‘Hubble’ flow case, we have shown that the solutions 
of the kinetic (Boltzmann) and hydrodynamic (Israel-Stewart) equations are essentially the 
same, in that they give the identical conserving energy momentum tensor. Most of the 
previous studies have found systematic differences between the solutions of the kinetic and 
hydrodynamic equations, the former is considered to be more fundamental. However, in the 
present case they agree exactly. 

We then considered the Boltzmann equation in the boost-invariant setup relevant to the 
hnal state of heavy-ion collisions. At asymptotically large times, the equation selects the 
Bjorken solution e ~ However, we have pointed out that in the conformally sym¬ 

metric case a novel solution of the form e = Cjr^ exists, and we precisely determined the 
normalization C which is not an arbitrary parameter. Very interestingly a solution with the 
same r-dependence was previously found as an exact solution of second order hydrodynam¬ 
ics. The existence of this solution is also suggested by fluid-gravity correspondence. But 
here our argument is not adequate, and more work is needed to firmly establish (or exclude) 
such a solution of the Einstein equation. 

We conclude by mentioning some possibilities to explore from our findings which might 
have potential applications in more realistic scenarios such as in high energy nuclear colli¬ 
sions. At early times after the collision between highly energetic nuclei, the produced plasma 
of quarks and gluons reaches high temperatures with almost vanishing chemical potential 
such that it is approximately conformal. This tiny plasma is far away from equilibrium due 
to the violent acceleration in all the spatial directions and thus, the 3D Hubble solution 
discussed in Sect. Hi] offers the possibility to model at this stage the system. In addition, the 
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3D Hubble solution can encode the early-time information about the flow configuration and 
energy density behavior ~ This scaling of the energy density at the beginning of the 
expansion might be connected with the unorthodox Bjorken flow solution whose energy den¬ 


sity decays faster ~ r than the standard Bjorken solution ~ r 2^. Therefore, within 


this scenario the lifetime of the fireball would be shorter since the unorthodox Bjorken flow 


predicts that the system reaches faster the freeze-out temperature. This possibility can be 


studied nowadays by performing a systematic hne-tuning analysis of the freeze-out tempera¬ 
ture, the initial time when the numerical hydrodynamical simulation starts to run and other 
parameters such as the transport coefficients. 
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Appendix A: Another derivation of (I40|) 

In order to be more confident with the unfamiliar number 13 = 1.27672 found in flTOj) . let 
us derive it from a slightly different perspective. Doing integration by parts in the right-hand 
side of fl38|) and keeping contributions only from t' = r, we get 




2(2 -/ 3 ) 2 io 


/3 (3 


I 


dx{x^ + 1)^ arctanx 
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where we changed variables as a; = — 1. The x-integral is convergent at x —)■ cxo 

when /9 > 3. We thus temporarily assume /? > 3 and send the upper limit to inhnity. Then 
the integral can be done exactly 


g 1, , d , (. P\ T(¥) , 3f2Q.i,l;|.3-|;l) \ 

r'-fl 2 { 2 -ff) 8 V r(|)r(3 - f) j 

We then analytically continue the result to /? < 3 and look for the zero of flA2p . 
jS ~ 1.27672 in perfect agreement with the previous method. 


(A2) 
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